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CHAPTER  1  INTRODUCTION  AND  CORRECTION  METHODOLOGIES 


1,1  Purpose 

Past  efforts  have  been  made  to  generate  a  suitable  wind  tunnel  bloekage  eorreetion 
methodology  using  wind  tunnel  wall  pressure  signatures.  Most  notable  among  these 
efforts  has  been  those  of  Glauert,  Masked,  Ashill  and  Keating,  Hensel,  Haekett- 
Wilsden,  and  Sorensen  and  Mikkelsen.  Most  of  the  methods  from  these  sourees  do 
not  make  use  of  measured  wall  pressures.  The  analyses  used  infer  what  tunnel  wall 
pressure  is  expeeted  to  be,  but  they  do  not  use  measured  values  exeept  for  Ashill  and 
Keating  and  Haekett-Wilsden.  It  was  the  purpose  of  this  thesis  to  eolleet  high  quality 
powered  propeller  data  for  a  wide  range  of  size  of  propeller  relative  to  the  wind 
tunnel  so  that  wall  eorreetion  methods  eould  be  evaluated  and,  perhaps,  extended.  To 
this  end,  powered  propeller  data  was  eolleeted  in  the  Glenn  L.  Martin  Wind  Tunnel  at 
the  University  of  Maryland  and  at  the  Naval  Aerodynamies  Test  Faeility  at  Patuxent 
River  Naval  Air  Station;  the  bloekage  of  the  first  was  negligible.  The  data  resulting 
from  this  test  was  taken  as  baseline  data.  The  bloekage  of  the  seeond  was  non- 
negligible  and  the  Haekett-Wilsden  and  Sorensen  and  Mikkelsen  eorreetion  methods 
and  Glauert  methods  were  employed,  separately,  to  examine  the  ealeulated  bloekage 
eorreetion  and  eompare  the  methodologies.  Critiques  and  eonelusions  are  drawn 
about  the  quality  of  this  eorreetion  methodology,  partieularly  as  it  applies  to  the 
powered  propeller  model. 
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1,2  Glauert  Correction 


Glauert  was  among  the  first  aero  dynamic  ists  to  explore  propeller  blockage 
corrections.  In  one  of  his  early  publications,  The  Elements  of  Aerofoil  and  Airscrew 
Theory,  Glauert  presents  the  derivation  of  his  correction.  This  correction  determines 
an  “equivalent  tree  airspeed,  V’,  corresponding  to  the  tunnel  datum  velocity  V,  at 
which  the  airscrew,  rotating  with  the  same  angular  velocity  as  in  the  tunnel,  would 
produce  the  same  thrust  and  torque.”'  This  equivalent  tree  airspeed  can  be  found  by 
applying  momentum  theory  to  the  case  of  an  airscrew  rotating  in  a  wind  tunnel.  In  a 
closed  jet,  under  positive  thrust  conditions,  the  propeller  slipstream  accelerates  and 
maintains  a  velocity  greater  than  the  wind  tunnel  datum  velocity,  V.  Due  to 
continuity,  the  same  volume  of  air  must  pass  just  ahead  of  the  propeller  as  passes  just 
aft  of  the  propeller.  Therefore,  it  follows  that  the  velocity  outside  the  slipstream  will 
be  less  than  V.  Due  to  its  lower  velocity,  the  static  pressure  outside  of  the  slipstream 
has  a  higher  static  pressure.  Using  Bernoulli  and  taking  into  account  a  pressure 
discontinuity  across  the  propeller  disc,  Glauert  shows  that  the  static  pressure  of  the  air 
within  the  slipstream  must  be  identical  to  that  outside  of  the  slipstream.  This  higher 
static  pressure  “reacts  back  to  the  propeller  so  that  it  develops  thrust  that  is  greater 
than  would  be  developed  in  an  unrestricted  flow  of  the  same  speed  with  the  same 
propeller  rotation  rate  and  blade  angle.  Or  it  can  also  be  said  the  thrust  developed 
would  be  equal  to  that  which  would  be  expected  at  a  lower  V’  in  tree  air.”  Glauert’s 
correction  follows: 

VA^’  =  [1  -  (T4ai  /  (2*(l+2*  T4)'^^))]''  (1) 


2 


Where  V’  is  the  equivalent  free  airspeed 

V  is  the  wind  tunnel  datum  veloeity 

14  =  T  /  (pAV")  (2) 

T  =  thrust,  lbs 
p  =  density,  Ibs/ft^ 

A  =  cross-sectional  area  of  propeller,  ft^ 

V  =  wind  tunnel  datum  velocity,  ft/s 
tti  =  A/C  (3) 

C  =  cross-sectional  area  of  wind  tunnel,  ft^ 

Glauert’s  results  are  shown  in  figure  1. 
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Speed  Correction  for  Thrusting  Propeller 


A/.A 


Figure  1.  Speed  correction  for  thrusting  propeller.  Glauert  Correction. 


It  can  be  seen  that  V’,  the  equivalent  free  airspeed,  is  less  than  V,  the  tunnel  datum 
speed,  for  any  positive  thrust  eonditions.  Negative  values  of  thrust  are  not  shown  as 
the  propeller  is  no  longer  in  a  positive  thrust  condition,  resulting  in  the  above 
slipstream  model  being  invalidated.  Glauert’s  correction  is  a  classic  and  simple 
eorreetion  that  is  used  in  modern-day  testing  to  eorreet  for  wind  tunnel  interferenee. 

1,3  Maskell  Correction 

Maskell  developed  a  theory  of  bloekage  eorreetion  for  bluff  bodies  in  elosed  wind 
tunnels  based  on  an  approximate  relation  deseribing  the  momentum  balanee  in  the 
flow  outside  the  wake  of  a  bluff  body  and  two  auxiliary  equations.  This  theory  was 
the  first  to  address  bloekage  eorreetions  for  non-streamline  flow,  ineluding  addressing 
flow  over  partially  stalled  shapes,  sueh  as  wings.  In  both  flow  oases,  streamline  and 
non-streamline  flow,  the  wake  aft  of  a  bluff  body  has  a  tendenoy  towards  axial 
symmetry.  Additionally,  through  experiments  oonduoted  at  the  Royal  Airoraft 
Establishment  (RAE),  it  was  shown  that  the  wall  oonstraint  oan  be  regarded  as 
equivalent  to  a  simple  inorease  in  velooity  of  the  undisturbed  stream.  Using 
oonservation  of  momentum  over  a  oontrol  surfaoe  whioh  is  bound  by  the  solid  walls 
of  the  wind  tunnel,  the  surfaoe  of  the  body,  and  the  oonstant-pressure  surfaoe 
bounding  the  effeotive  wake  and  two  planes  normal  to  the  undisturbed  velooity 
veotor — one  lying  upstream  of  the  body  and  the  other  looated  where  the  oross- 
seotional  dimensions  of  the  bubble  are  greatest— Maskell  developed  a  relation  for  the 
drag  ooefficient  based  on  base-pressure  measurements  of  the  model  in  the  tunnel^. 
This  relationship  is  affirmed  by  further  RAE  experiments.  Maskell  then  uses 
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relationships  describing  the  distortion  of  the  wake  due  to  wall  constraint,  combined 
with  the  drag  coefficient  relationship  to  generate  a  blockage  correction  expression. 
This  expression  is 

Aq/q  =  e  (CdS/C)  (4) 

Where  e  is  the  blockage  factor 
Cd  is  the  drag  coefficient 
S  is  the  reference  area  of  the  model 
C  is  the  cross-sectional  area  of  the  wind  tunnel 
Aq  =  qc  -  q  is  the  effective  increase  in  dynamic  pressure  of  the 
undisturbed  stream  due  to  constraint 

The  blockage  factor  is  dependent  on  the  base-pressure  coefficient 


e=l/(Ke"-l)  (5) 

Where  Kc  =  1-Cpbc  and  is  found  iteratively 

Cpbc  is  the  corrected  base-pressure  coefficient 

“The  factor  e  is  shown  to  range  between  a  value  of  a  little  greater  than  5/2  for  axi- 
symmetric  flow  to  a  little  less  than  unity  for  two-dimensional  flow.  But  the  variation 
from  5/2  is  found  to  be  small  for  aspect  ratios  in  the  range  of  1  to  10.^”  To  show  the 
validity  of  the  theory,  data  was  collected  for  the  two  extremes  for  which  this  theory 
purportedly  is  valid:  axi-symmetric  flow  and  two-dimensional  flow.  The  results  show 
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that  “in  the  interval  A=(l,10)  the  bloekage  faetor  lies  roughly  in  the  range  e  =  5/2  ± 

V4.  And  the  eonstant  value  e=5/2  leads  to  errors  of  only  ±  0.1  Aq  at  the  extreme  points 
of  the  range.”  “Therefore,  the  theory  holds  for  almost  all  two-dimensional  bluff-body 
flows,  and  for  the  wide  range  of  three-dimensional  flows  for  whieh  the  wake  is 
elosely  axi-symmetrie  at  the  downstream  plane.”  However,  it  must  be  noted  that 
there  is  an  important  possible  exeeption  to  this  rule.  An  implied  assumption  in  the 
theory  is  that  the  origin  of  the  wake  is  independent  of  eonstraint.  And  so  it  may  be 
neeessary  to  exelude  well-rounded  bluff  bodies  for  whieh  a  small  ehange  in  pressure 
distribution  might  lead  to  a  signifieant  movement  of  the  separation  front.  This  also 
limits  the  eonfidence  in  correction  of  partially-separated  flows  such  as  those  over 
stalled  wings. 

MaskelTs  method  is  not  designed  for,  or  easily  adaptable  to,  application  of  powered 
propellers  and  therefore  this  methodology  was  not  used  for  blockage  corrections  in 
this  powered  propeller  test.  Its  description  is  included  here  because  it  uses  a 
measured  value  of  drag  to  infer  an  addition  to  blockage  of  the  solid  object  in  the  test 
section.  This  is  somewhat  analogous  to  the  Glauert  propeller  blockage  approach. 

1,4  Ashill  and  Keating  Correction 

Ashill  and  Keating,  of  the  Royal  Aircraft  Establishment,  build  off  of  MaskelTs  work 
in  development  of  their  correction  theory.  As  Maskell  did,  Ashill  and  Keating  create 
a  method  of  correction  which  does  not  require  modeling  of  flows  in  the  tunnel.  This 
eliminates  the  difficulty  of  modeling  separated  and  complicated  flows.  Instead,  they 
use  a  ‘two-component’  method  in  which  the  knowledge  of  two  components  of  flow 
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velocity  at  the  outer  boundary  must  be  known.  These  components  are  that  normal  to 
the  boundary  and  that  in  the  stream  wise  direction.  For  low  speed  wind  tunnels  that 
have  solid  walls  and  thin  wall  boundary  layers,  the  outer  boundary  of  the  ‘two- 
component’  method  may  be  taken  to  coincide  with  the  tunnel  walls  and  the  normal 
component  assumed  to  be  zero.  Hence  only  the  stream  wise  component  is  required 
for  tunnels  of  this  type,  and  this  component  can  be  inferred  from  static -pressure 
measurements  at  the  walls"*. 

Using  the  two-component  method  as  described  above,  Ashill  and  Keating  use 
equations  for  the  perturbation  potential  in  free  air  flow  and  that  due  to  wall  constraint. 
They  show  that  the  difference  in  perturbation  potential  between  free  air  and  the 
constrained  condition  is  entirely  due  to  wall  interference.  They  use  Bernoulli’s 
equation  to  calculate  the  stream  wise  velocity  increment  from  the  corresponding 
increment  in  static  pressure.  This  increment  can  be  used  to  asses  the  difference  in 
perturbation  potential,  and  thus  wall  interference. 

To  test  this  theory,  Ashill  and  Keating  test  two  types  of  plates  in  the  RAE  tunnel;  one 
of  aspect  ratio  of  one  and  the  other  with  aspect  ratio  of  three.  They  also  test  a 
combat-aircraft  half-model.  In  each  case,  wall  static -pressure  taps  were  used  to  asses 
changes  in  pressure  along  the  wall.  Base-pressure  coefficients  are  also  measured. 

This  data  allows  for  calculation  of  corrected  base-pressure  coefficients  and  blockage 
effects"*.  The  comparisons  of  Ashill  and  Keating’s  corrected  base-pressure 
coefficients  against  MaskelTs  show  good  agreement  with  MaskelTs  theory  for  the  flat 
plates.  However,  for  the  combat-aircraft  half-model  the  blockage  corrections  are 
significantly  smaller  than  those  given  by  MaskelTs  method. 
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As  with  Maskell’s  method,  Ashill  and  Keating’s  correction  methodology  was  not 
developed  with  powered  propeller  testing  and  blockage  effects  in  mind  and  its 
application  to  propeller  corrections  has  not  been  reported  in  the  literature.  It  could  be 
adapted  to  propeller  corrections  using  an  approach  analogous  to  the  that  used  in  this 
study  with  the  Hackett-Wilsden  method  as  discussed  in  later  sections. 

1.5  Hensel’s  Area-Ratio  Method 

Hensel’s  1951  technical  note  reports  a  method  for  calculating  the  “ratios  of  the 
velocity  increments  at  test  bodies  to  those  at  the  tunnel  walls  caused  by  the  solid 
blocking  of  these  bodies  within  the  walls  of  a  closed  rectangular  wind  tunnel.”^ 
Hensel  builds  strongly  off  the  methods  in  Thom^  in  which  point  and  line  sources  are 
used  to  model  bodies  of  revolution,  finite  straight  wings,  and  finite  swept  wings. 

1.6  Hackett-Wilsden  Method 

In  1979,  Hackett,  Wilsden,  and  Lilley  released  a  NASA  report^  describing  an 
expansion  on  Hensel’ s  correction  methodology.  Instead  of  using  doublets  to  create 
an  “equivalent  body”  outline,  they  used  sources  and  sinks  to  do  so.  In  their 
experiments,  they  used  models  that  produced  varying  blockage  ratios;  some  models 
were  small  enough  such  that  they  were  in  a  “free  air”  state.  The  largest  models  were 
up  to  10%  of  the  total  tunnel  cross  section.  During  these  experiments,  they  measured 
wall  pressures  along  the  center  line  of  the  tunnel  walls  or  the  roof,  and  on  the  floor 
when  necessary.  In  cases  where  the  body  created  lift,  the  mean  of  the  roof  and  floor 
pressures  was  used.  These  pressures  were  then  “used  to  determine  source  or  sink 
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strengths,  spans,  and  locations  on  the  tunnel  center  line  which  define  a  body  outline 
which  ‘is  equivalent’  to  the  test  model  and  its  wake7” 

The  Hackett-Wilsden  method  employs  an  iterative  scheme  to  resolve  the  pressure 
signatures  into  symmetric  and  anti-symmetric  parts.  The  symmetric  part  represents 
the  solid  blockage  while  the  anti-symmetric  part  reflects  the  wake  blockage.  Once 
these  parts  are  resolved  and  the  necessary  parameters  are  deduced  Irom  the  parts,  a 
set  of  equations  and  a  chart  look-up  technique  is  employed  to  obtain  singularity 
strengths,  spans,  and  locations.  Using  these  parameters,  a  velocity  increment  Au/Uoo 
is  calculated.  This  is  the  velocity  increment  calculated  Ifom  the  singularity  strengths, 
spans,  and  locations  which  represent  an  “equivalent  body”  outline;  therefore,  this  is 
also  the  velocity  increment  caused  by  the  model  in  the  tunnel. 

Using  the  calculated  velocity  increment,  the  Hackett-Wilsden  method  then  employs 
one  of  two  correction  methods.  For  models  with  predominately  attached  flow,  a 
global  correction  can  be  used.  To  do  this,  a  source-panel  model  of  the  body  may  be 
generated  and  the  three  components  of  interference  velocity  at  each  panel  center  is 
calculated.  These  are  then  input  to  a  source  panel  program  together  with  the  panel 
geometry  and  mainstream  velocity  to  obtain  corrected  pressure  coefficients. 

For  models  that  do  not  have  predominately  attached  flow,  a  point-by-point  correction 
equation  is  used  to  obtain  corrected  pressure  coefficients.  It  follows: 


Cpe(x)  =  [(Cpu(x)  -  1)  /  (1  +  (  Au(x)/Uoo)'  )]  +  1  (6) 
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To  prove  the  validity  of  the  point-by-point  eorreetion,  in-tunnel  pressure  eoeffieients 
taken  with  a  cheekout  sphere  installed  were  eorreeted  using  the  above  equation  and 
eompared  with  a  eorresponding  free-air  ealeulation.  The  resulting  error  in  Cp  was 
0.013,  at  most.  Along  the  equator  of  the  sphere,  the  error  was  less  than  0.005  at  all 
positions.  Though  the  bloekage  ratio  for  this  eheekout  sphere  is  not  explieitly  stated, 
results  do  state  that  bearing  in  mind  the  large  size  of  the  model  used  for  eheekout,  it 
may  be  eoneluded  that  the  above  equation  eorreets  measured  pressures  suffieiently 
accurately  for  most  practical  purposes. 

For  instance,  for  models  occupying  up  to  13.7%  of  the  tunnel  cross  section,  the 
Hackett-Wilsden  methodology  produced  results  comparable  with,  and  no  less  credible 
than,  results  from  similar  test  in  a  large  wind  tunnel.  This  allows  for  more  confidence 
in  conducting  blockage  corrections  for  relatively  large  models.  Additionally,  for 
models  with  predominately  attached  flow,  a  global  correction  may  be  used 
successfully. 

1,7  Sorensen  and  Mikkelsen  Method 

In  2002,  Sorensen  and  Mikkelsen  presented  a  method  to  model  wind  tunnel  blockage 
with  particular  attention  to  the  testing  of  wind  turbines.  This  method  is  based  on  the 
axial  momentum  theory  and  employs  a  series  of  five  equations  to  obtain  a  closed- 

o 

form  solution  of  the  equations  of  the  axial  momentum  theory  .  It  incorporates 
Hackett’s  method  of  using  source  and  sink  singularities  to  model  the  contraction  or 
expansion  of  the  wake  of  the  propeller.  Using  the  strength  of  the  specific  singularity 
which  accurately  models  this  expansion  or  contraction,  Sorensen  and  Mikkelsen 
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introduce  an  explicit  formula  for  calculating  an  intermediate  value,  u,  such  that  the 
velocity  through  the  propeller  disc  can  be  calculated  directly  and  all  five  equations 
derived  Ifom  the  axial  momentum  theory  can  be  solved  simultaneously.  This 
intermediate  value  is  also  used  to  obtain  the  correction  synonymous  to  Glauert’s 
correction  formula.  This  explicit  equation,  as  well  as  the  five  momentum  equations, 
and  the  Sorensen  and  Mikkelsen’s  correction  formula  are  shown  below: 

Intermediate  Variable  Equation 

u  =  (a  (Pa^  -  1))  /  (Pa  (3a  -  2)  -  2a  +  1) 

Where  P  =  S/C 

S  =  propeller  cross  sectional  area,  ft^ 

C  =  tunnel  cross  sectional  area,  ft^ 

a  =  strength  of  singularity,  Si/S 

Si  =  cross  sectional  area  of  slip  stream 
far  downstream,  ft^ 

u  =  u/Vo 

u  =  velocity  of  air  through  propeller,  ft/sec 
Vo  =  wind  tunnel  datum  velocity,  ft/sec 

Five  Momentum  Equations 

ui  =  ou  (10) 

Where  ui  =  Ui/Vo 

ui  =  slip  stream  airspeed  far  downstream,  ft/sec 
U2(l-pa)=  1 -pu  (11) 


(V) 

(8) 
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Where  U2  =  U2W0 


U2  =  airspeed  outside  of  slip  stream 
far  downstream,  ft/see 

Ct  =  Ui^-U2^  (12) 

Where  Cj  =  coeffieient  of  thrust 

Cp=l-U2'  (13) 

Where  Cp  =  eoeffieient  of  power 

-pCi  -  Cp  =  2uipa(ui  -  1)  -  2u2(l-pa)(l  -  U2)  (14) 

Correetion  Formula 

VWo  =  u-(y4*CT/u)  (15) 

Where  V’  is  the  equivalent  free  airspeed 

1,8  Summary 

Chapter  2  of  this  doeument  deseribes  the  experimental  set  up  of  the  test,  deseription 
of  test  facilities,  and  the  steps  taken  to  reduce  the  data  from  the  experiment.  Chapter 
3  gives  detailed  analysis  steps  and  methodology  that  was  used  to  first  mathematically 
simulate  the  physical  phenomena  in  the  wind  tunnel  using  the  Hackett-Wilsden 
approach  and  to  calculate  the  Hackett-Wilsden  and  Sorensen  and  Mikkelsen 
corrections.  The  Glauert  correction  is  also  presented.  Pressure  predictions  using 
Glauert  and  Sorensen  and  Mikkelsen  equations  are  explained.  Discussion  on  the 
corrections  and  the  pressure  predictions  is  given.  Finally,  Chapter  4  offers 
conclusions  and  Chapter  5  proposes  recommendations  for  future  work. 
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It  should  be  noted  that  for  eaeh  propeller,  nine  sueeessful  eonditions  were  desired  for 
exeeution  of  analysis.  These  eonditions,  eorresponding  to  three  wind  tunnel  speeds, 
eaeh  at  three  RPM  settings,  would  be  simulated  using  sourees  and  sinks.  With  nine 
runs  for  eaeh  of  four  propellers,  a  total  of  thirty-six  potential  runs  could  be  analyzed 
and  mathematically  simulated.  However,  in  some  cases,  high  speed/high  RPM  runs 
were  unable  to  be  completed  successfully;  the  motor  overheated  prior  to  the 
successful  completion  of  the  run.  As  many  consistent  conditions  as  possible  were 
identified  and  analyzed. 

In  the  body  of  this  document,  experimental  data,  calculations,  corrections,  and 
predictions  are  shown  for  only  the  I8”x8”  propeller  runs.  Other  propeller  data, 
calculations,  etc,  are  not  shown  unless  specific  attention  is  required  in  analysis  or 
explanation.  All  data  for  the  I4”x8”,  22”x8”,  and  24”x8”  propellers  can  be  found  in 
Appendix  B. 
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CHAPTER  2  DESCRIPTION  OF  TEST 


2.1  Description  of  Test  Articles 

2.1.1  Propellers 

The  Scimitar  Series  propellers  used  were  acquired  from  Windsor  Propeller  Company. 
The  four  propellers  had  diameters  of  14”,  18”,  22”,  and  24”  and  identical  airfoil 
shapes.  All  propellers  are  fabricated  from  hard  wood  maple.  A  drill  guide  balance 
system,  purchased  from  Windsor  Propeller  Company,  was  used  to  balance  the 
propellers.  The  balance  stand  allowed  tip-to-tip  and  chord-wise  balancing,  while  the 
drill  guides  and  bushings  provided  a  guide  though  which  to  drill  the  centered 
mounting  hole,  using  the  included  3/8”  drill  bit.  The  figures  in  Appendix  A  show  the 
experimental  set  up  and  articles. 

2. 1 .2  Motor  and  Electronic  Speed  Controller 

An  AXl  5330-18  outrunner  brushless  motor  was  used  to  drive  the  propellers.  The 
motor  is  2.5”  in  diameter  and  2”  in  length.  The  maximum  power  output  is 
approximately  3.5  horsepower  (2.6  kW).  The  propeller  shaft  was  12  mm  in  diameter 
(0.47”)  when  purchased,  but  was  turned  down  to  3/8”  to  mount  the  propellers.  The 
motor  is  controlled  by  a  Jeti  Advance  PLUS  90  Opto  Control  brushless  electronic 
speed  controller  (ESC).  The  remote  control  receiver  was  connected  to  the  ESC, 
allowing  the  ESC  to  be  controlled  by  the  remote  control  exterior  to  the  tunnel. 
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2.1.3  Torquemeter 

The  torque  cell  used  during  testing  was  a  5330  Hollow  Flange  Reaction  Torque  Cell, 
manufactured  by  Interface.  Its  maximum  capacity  is  60  in- lb.  The  diameter  of  the 
torque  cell  is  2.5”  in  diameter  and  2.125”  in  length. 

2.1.4  Power  Supply 

The  power  supply  was  rented  from  Elgar  Electronics.  It  requires  200/208  volt,  three- 
phase  input  and  is  capable  of  providing  4  kW  of  power.  The  maximum  voltage  was 
set  for  42.3  volts  and  the  amperage  was  varied  from  0  to  50  amps  based  on  demand 
from  the  ESC,  which  was  controlled  by  the  remote  control.  This  voltage  and 
amperage  limitation  was  introduced  so  as  not  to  exceed  the  motor  limitation  of  2kW 
of  power  input. 

2.1.5  Remote  Control 

A  typical  2-channel  RC  car  or  plane  remote  control  was  used.  It  was  bought  from 
Tower  Hobbies  and  included  the  necessary  signal  receiver.  Modifications  were  made 
to  channel  2,  throttle.  The  throttle  stick  was  replaced  with  a  trim  pot  so  that  a  power 
setting  could  be  set  and  the  knob  released  without  the  setting  changing. 

2.1.6  Stroboscope 

A  1531 A  type  stroboscope,  manufactured  by  General  Radio  and  owned  by  the  Glenn 
E.  Martin  Wind  Tunnel  (GEMWT)  at  the  University  of  Maryland  was  used  to  assess 
the  RPM  of  the  propellers  during  testing. 
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2.1.7  Pressure  Board 


A  pressure  board,  owned  by  GLMWT,  was  mounted  on  the  ceiling  of  each  tunnel  to 
record  the  pressure  signature  during  running.  In  the  GLMWT,  two  pressure  boards 
were  used.  These  boards,  together,  are  12’  long  and  have  24  pressure  ports  spaced 
every  6”  along  the  center  of  the  board.  It  was  mounted  with  bolts  and  sealed  on  the 
sides  with  aluminum  tape.  The  pressure  tubes  were  taped  to  the  side  of  the  tunnel  and 
run  through  the  floor  to  the  sensing  equipment. 

One  of  these  two  boards  was  transported  for  use  at  NATF.  The  board  used  is  7’  long 
and  has  14  pressure  ports  spaced  every  6”  along  the  center  of  the  board.  It  was  also 
mounted  on  the  center  of  the  wind  tunnel  test  section  ceiling  and  sealed  with 
aluminum  tape.  Pressure  tubes  were  routed  out  of  the  ceiling  of  the  wind  tunnel. 

2.2  Experimental  Setup 

2.2.1  Support  Structure 

The  motor  was  attached  to  a  bearing  and  housing  which  was  fabricated  to  prevent 
cantilevered  loads  on  the  torque  cell.  This  housing  was  attached  to  a  horizontal 
support,  which  was  then  mounted  to  a  vertical  circular  steel  support.  For  use  in  the 
GLMWT,  this  support  was  approximately  4’  tall  which  allowed  for  the  propeller  to  be 
in  the  center  of  the  tunnel.  An  aerodynamic  fairing  was  used  to  cover  the  bottom  two 
feet  of  this  support  to  aid  in  drag  reduction  and  shield  the  more  bulky  remote  control 
receiver.  This  support  was  mounted  to  the  tunnel  balance. 

The  three  power  wires  protruding  from  the  motor  were  soldered  to  the  three  wires 
coming  from  the  ESC.  The  ground  lines  from  each  were  soldered  together,  and  the 
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remaining  wires  were  attached  such  that  the  motor  turned  in  a  counter  clockwise 
direction  (from  the  front).  The  wires  protruding  from  the  opposite  side  of  the  ESC 
were  soldered  to  longer  lengths  of  wire  fed  through  the  floor  of  the  tunnel.  These 
lines  were  attached  to  the  output  of  the  power  source.  The  remote  control  receiver 
was  plugged  into  the  ribbon  wire  coming  from  the  ESC. 

A  terminal  block  was  mounted  on  the  back  of  the  housing.  This  was  used  to  connect 
the  torque  cell  output  wires  to  the  wires  fed  into  the  instrumentation,  which  were  also 
fed  through  the  floor  of  the  tunnel. 

The  full  support  structure,  as  installed  in  the  GEMWT  and  Naval  Aerodynamic  Test 
Eacility  (NATE)  are  show  in  figures  in  Appendix  B. 

2,3  Description  of  Wind  Tunnels 

2.3.1  Physical  Characteristics 

The  Glen  L.  Martin  Wind  Tunnel  at  the  University  of  Maryland  and  the  Naval 
Aerodynamic  Test  Eacility  Subsonic  Wind  Tunnel  at  Patuxent  River  Naval  Air 
Station,  MD,  were  used  for  testing.  The  GEMWT  is  a  closed-return  tunnel  and  has  a 
7.75H  X  1 1.04W  foot  test  section  that  is  13.25  feet  long.  It  has  corner  fillets  and  a 
turbulence  intensity  factor  of  0.21%.  The  test  section  area  is  85.04  sq  feet.  The  speed 
range  is  from  3  ft/s  to  337  ft/s. 

The  Naval  Aerodynamic  Test  Eacility  Subsonic  Wind  Tunnel  at  Patuxent  River  Naval 
Air  Station  has  a  4  x  4  foot  test  section  that  is  8  feet  long.  It  has  sharp  corners  and  a 
turbulence  intensity  factor  of  0.6%.  The  test  section  area  is  16  sq  feet.  The  speed 
range  is  from  40  ft/s  to  200  ft/s. 


18 


For  data  reduction  purposes,  dimensions  of  the  NATF  tunnel  were  normalized  in 
several  ways.  For  modeling  of  solid  blockage,  calculation  of  velocity  increment  at 
the  tunnel  ceiling,  and  method  of  images  calculations,  wind  tunnel  dimensions  were 
normalized  by  the  test  section  span.  The  position  of  x/B  =  0  was  assigned  to  the 
center  of  the  model  length,  slightly  upstream  of  the  mid-length  position  in  the  tunnel. 
Using  this  position,  the  beginning  of  the  wind  tunnel  test  section  occurs  at  x/B  =  -0.85 
and  ends  at  x/B  =  1.148.  The  height  of  the  test  section  is  z/B  =  -0.5  to  z/B  =  0.5  and 
the  width  of  the  test  section  is  y/B  =  -0.5  to  y/B  =  0.5. 

To  provide  comparison  of  measured  pressure  data,  experimental  results  were 
normalized  by  propeller  diameter,  D.  The  x/D  =  0  was  assigned  to  the  location  of  the 
propeller  in  this  case. 


2.3.2  Data  Acquisition  Hardware 

GLMWT  uses  a  six  component  yoke  type  balance.  Pressure  measurements  are 
obtained  using  an  Electronically  Scanned  Pressure  system  made  by  Pressure  System 
Inc.,  coupled  with  pressure  boards,  as  described  previously.  Up  to  144  static  pressure 
measurements  may  be  obtained.  Twenty- four  measurements  were  obtained  for  this 
test. 

NATF  uses  a  three  component  platform  type  balance  to  measure  side  force,  yaw 
force,  and  axial  force.  Pressure  measurements  are  obtained  using  a  system  identical 
to  that  at  GLMWT,  the  PSI  8400  Electronically  Scanned  Pressure  system  made  by 
Pressure  System,  Inc.,  coupled  with  the  pressure  board  mounted  at  NATE.  The 
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pressure  system  is  capable  of  taking  up  to  256  pressure  port  measurements.  Fourteen 
pressure  measurements  were  obtained  for  this  test. 

NATF  has  a  Hitachi  HV-D5W  digitahanalog  broadcast  quality  camera  that  can  be 
placed  for  optimum  visualization.  A  Panasonic  analog  camera  is  also  mounted 
looking  into  the  test  section  from  the  top. 

2.3.3  Data  Acquisition  Software 

GLMWT  uses  a  Hewlett  Packard  743  VXI  bus  workstation.  Customized  data 
acquisition  systems  and  graphical  user  interfaces  may  be  designed  using  icon-driven 
software  VEETEST. 

The  data  acquisition  software  and  systems  at  NATE  are  dependent  on  the  data 
acquisition  instrumentation  used.  Eor  this  experiment,  a  486  PC  Autonet  system  was 
used. 


2,4  Test  Execution 

Three  RPM  settings  were  used  for  testing  each  propeller.  These  RPM  settings 
corresponded  to  set  current  on  the  power  supply.  The  power  supply  was  set  for  a 
maximum  voltage  of  42.3  at  all  times  due  to  the  model  motor  limitations.  Using  the 
remote  control,  the  amperage  output  was  selected  and  set,  a  complete  wind  tunnel 
velocity  sweep  was  conducted  with  stabilized  points  for  data  collection,  and  then  the 
amperage  was  increased  or  decreased  for  a  different  RPM  setting.  Three  amperage 
settings  were  used  to  obtain  three  RPM  settings:  low,  medium,  and  high.  Because  the 
power  required  to  rotate  the  propeller  was  greater  at  low  wind  tunnel  speeds  than  high 
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speeds,  eurrent  settings,  partieularly  the  high  eurrent  setting,  were  set  at  a  low  tunnel 
speed.  This  guaranteed  that  the  power  required  by  the  motor  did  not  exceed  the  2kW 
from  the  power  supply  as  the  tunnel  speed  varied  and  allowed  for  a  full  range  of 
tunnel  speeds  to  be  tested  at  a  consistent  power  setting. 

The  trim  pot  on  the  modified  remote  control  proved  to  be  very  sensitive  and  it  was 
difficult  to  set  power  supply  current  in  a  repetitive  and  efficient  fashion. 

Additionally,  the  actual  power  required  by  the  propeller,  for  a  given  setting  of  the 
remote  trim  pot,  increased  as  wind  tunnel  speed  decreased.  Therefore,  because  of  this 
sensitivity,  the  change  in  required  power,  and  for  efficiency,  tests  were  conducted  by 
setting  a  current  setting  with  the  remote  control  and  varying  wind  tunnel  speed.  First, 
the  current  was  set  at  the  lowest  setting  possible  with  the  wind  tunnel  velocity  at  a 
high  speed.  This  is  the  minimum  power  required  by  each  propeller  in  the  sequence  of 
testing  and  resulted  in  the  low  RPM  setting.  The  low  current  setting  was  maintained 
as  the  tunnel  velocity  was  decreased.  While  the  setting  on  the  remote  control  was 
never  changed,  the  power  required  by  the  motor  increased  and  the  curent  drawn  from 
the  power  supply  increased.  The  resulting  current  required  at  the  low  wind  tunnel 
speed  then  defined  the  “low  power  setting”  boundary.  By  knowing  that  the  maximum 
power  output  of  the  power  supply  was  limited  to  2kW  and  the  maximum  voltage 
setting  was  42.3  volts,  limited  by  the  motor,  the  maximum  current  available  could  be 
calculated  as  50  amps.  With  this  knowledge,  an  intermediate  current  setting  was  then 
deduced.  The  low  wind  tunnel  speed  was  maintained  as  the  intermediate  current  was 
applied,  and  thus  the  intermediate  RPM  setting  was  obtained.  Time  was  given  for  the 
propeller  to  attain  a  constant  RPM  and  then  data  was  taken  at  the  intermediate  RPM 
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setting  as  the  tunnel  was  increased  in  speed.  Finally,  because  the  power  required  was 
greatest  at  low  wind  tunnel  speeds,  the  tunnel  speed  was  brought  back  down  to  a  low 
speed  and  the  maximum  current  of  50amps  was  set.  The  tunnel  was  again  increased 
in  speed  as  data  was  taken  for  the  high  RPM  setting. 

To  define  the  highest  wind  tunnel  speed  at  which  data  was  taken  for  each  of  the  RPM 
settings,  the  tunnel  speed  was  varied  Ifom  the  minimum  speed  to  a  speed  such  that  the 
propeller  no  longer  produced  thrust.  This  was  determined  by  noting  the  measured 
drag  at  each  test  point,  subtracting  the  appropriate  drag  tare  value,  and  determining  if 
zero  thrust  condition  had  been  reached.  The  tunnel  speed  was  increased  and  test 
points  taken  until  this  occurred. 

During  each  test  point  measurements  of  drag  and  torque  was  recorded.  The  pressure 
at  each  tap  along  the  pressure  board  was  also  recorded. 

Wind  tunnel  testing  was  first  conducted  at  GLMWT.  Testing  was  begun  with  a  drag 
tare  run.  The  18”  x  8”  propeller  was  run  first,  and  was  followed  by  the  22”  x  8”  and 
24”  X  8”  propellers.  Next  the  14”  x  8”  propeller  was  tested  to  complete  the  8”  pitch 
propeller  group.  The  18”  x  8”  propeller  was  re-tested  to  provide  for  examination  of 
repeatability. 

At  NATF,  the  first  run  was  also  a  drag  tare  run.  Again,  the  18”  x  8”  propeller  was 
tested  first.  It  was  tested  as  a  “check  run”  to  compare  against  the  GLMWT  running. 
When  the  test  was  completed  successfully  and  the  data  examined,  it  was  determined 
that  both  the  installation  and  data  collected  was  sound.  The  experiment  proceeded 
with  the  testing  of  the  14”  x  8”,  22”  x  8”  and  24”  x  8”  propellers.  For  the  24”  x  8” 
propeller,  a  static  thrust  run  was  also  completed.  That  is,  the  tunnel  was  left  turned 
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off  and  data  was  collected  as  the  current  setting,  and  thus  propeller  RPM,  was  varied. 
This  completed  the  testing  of  the  8”  pitch  propeller  group. 

Finally,  for  repeatability,  the  18”  x  8”  propeller  was  re-tested.  A  static  thrust  run  and 
windmill  run  was  included.  The  final  run  was  the  drag  tare  run  with  the  motor 
spinning.  This  concluded  the  experiment. 


2,5  Data  Reduction 

2.5.1  Wind  Off  Zero 

Wind  off  zero  (WOZ)  points  were  taken  before  each  run  at  GLMWT.  The  data 
acquisition  software  at  GLMWT  corrected  for  the  WOZ  points  automatically.  The 
data  that  was  obtained  by  the  experimenter  included  this  correction. 

Wind  off  zero  points  were  taken  before  and  after  eaeh  run  at  NATF.  Drag  values 
varied  insignificantly  between  the  WOZ  points  before  and  after  the  respective  run. 
To  eorrect  for  the  WOZ  drag  values,  the  first  and  last  points  were  averaged  and 
subsequently  subtracted  Ifom  eaeh  drag  value  eollected  during  the  run. 

In  some  cases,  there  was  a  significant  difference  in  torque  values  between  the  WOZ 
points  before  and  after  the  each  run.  A  “step-methodology”  was  used.  That  is,  an 
equal  step  was  added  or  subtracted  to  the  first  WOZ  torque  value  at  each  point  such 
that  the  last  WOZ  torque  value  was  reaehed  by  the  end  of  the  run.  This  value  was 
subsequently  subtracted  Ifom  each  torque  value  collected  during  the  run.  This 
provided  for  an  equal  change  in  WOZ  torque  value  as  the  run  progressed. 

2.5.2  Drag  Tare 
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A  drag  tares  was  taken  at  both  GLMWT  and  NATF.  At  both  tunnels,  these  were 
taken  without  the  motor  running  at  a  range  of  tunnel  dynamie  pressures.  A  drag  tare 
run  was  also  taken  at  NATF  with  the  motor  running  to  determine  if  there  was  an 
appreeiable  difference  in  drag  due  to  the  spin  of  the  motor.  As  shown  in  the  first 
figure  in  Appendix  B,  there  was  no  appreciable  difference  in  these  drag  tares.  The 
drag  tares  were  corrected  for  WOZ  points  as  described  above. 


2.5.3  Pressure  Signatures 

The  pressure  signature  measurements  required  reduction  in  order  to  use  them  with  the 
wall  interference  correction  methodology.  The  raw  data  received  ifom  NATF  was 
the  difference  between  the  pressure  measured  at  the  pressure  board  tap  and  the 
reference  pressure.  The  reference  pressure  was  the  barometric  pressure  of  the  day. 
This  data  is  divided  by  the  dynamic  pressure  in  the  tunnel  to  obtain  the  coefficients  of 
pressure  at  the  ceiling  along  the  length  of  the  tunnel. 

For  comparison  to  the  simulated  model  and  wake,  the  pressure  signatures  were 
converted  to  velocity  increments.  To  do  so,  the  pressure  signature  was  measured  at  a 
range  of  dynamic  pressures  without  the  model  present  in  the  wind  tunnel.  These  were 
converted  to  coefficients  of  pressure,  as  stated  above.  Then,  the  difference  between 
the  coefficients  of  pressure  with  the  model  and  without  the  model,  at  approximately 
the  same  dynamic  pressure,  was  calculated.  Finally,  this  was  turned  into  a  velocity 
increment  using  Wilsden/Hacketf  s  equation^  as  below: 
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Au/Uoo  =  (l-ACp)'^^-l 


(16) 


This  quantity  is  used  later  for  eomparison  purposes. 


2,6  Data  Uncertainty 

Uneertainty  exists  in  all  instrumentation  used  in  the  test  set  up.  At  GLMWT,  the 
uncertainty  in  the  pressure  data  is  ±0.0015  psi  and  the  uncertainty  in  the  drag  data, 
measured  by  the  yoke  balance,  is  ±0.015  pounds.  More  information  on  the  GLMWT 
and  its  data  acquisition  system  can  be  found  at  http://windvane.umd.edu/index.php. 

At  NATF,  the  uncertainty  in  the  pressure  data  is  ±0.1%  of  the  full  scale  value  and  the 
uncertainty  drag  data  is  ±0.818%  of  the  full  scale  value,  corresponding  to  ±0.49 
pounds.  The  uncertainty  in  the  tunnel  velocity  is  0.75%  of  the  full  scale  value, 
corresponding  to  1.5  ft/s.  More  information  on  the  NATF  tunnel  and  its  data 
acquisition  equipment  can  be  found  at  http  ://sata.tamu.  edu/ members/tunnels/ 1 3  8  .html 
and  in  reference  9  and  10. 

The  torquemeter  mounted  on  the  model  also  has  uncertainty.  According  to  the 
manufacturer’s  data  sheet,  the  uncertainty  in  torque  measurements  is  0.25%  full  scale, 
corresponding  to  0.15  lb-in. 
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CHAPTER  3  MATHEMATICAL  SIMULATION  OF  PHYSICAL 
PHENOMENA 


3.1  Model  and  Wake  Simulation 

3.1.1  Solid  Body  Simulation 

The  volume  of  the  model,  as  mounted  at  NAIF,  was  simulated  by  a  three- 
dimensional  source  and  sink  pair.  The  pair  was  placed  on  the  centerline  of  the  tunnel, 
as  the  model  was.  This  results  in  the  position  of  any  source  or  sink  [xq,  yo,  Zq]  equal 
to  [xo,  0,  0],  where  Xq  gives  the  position  of  the  singularity  along  the  test  length  of  the 
test  section.  The  source  and  sink  pair  is  superimposed  with  a  one-dimensional 
uniform  velocity,  which  simulates  the  free  stream  flow  inside  the  wind  tunnel. 
Derived  from  Katz  and  Plotkin",  the  velocity  potential  for  a  free  stream  flow 
superimposed  with  a  source  of  strength  Oss  at  position  [x,y,z]  =  [-Xo,0,0]  and  a  sink  of 
strength  Osk  at  position  [x,y,z]  =  [xo,0,0]  is 

0(x,y,z)  =  UooX-  Oss/  {(47i)[(x+Xo)^+(y)^+(zo)]'^^} 

+  ask/  {(47i)[(x-Xo)^+(yf+(z)]''2}  (17) 


The  stream  function,  with  cylindrical  coordinates  of  the  source  and  sink  of  [x,r]  =  [- 
Xo,0]  and  [x,r]  =  [xo,0],  respectively,  is 

¥(x,r)  =  Uoor  -  Gss/  {2[(x+Xo)2+(r)2]''2} 

ask/  {2[(x-Xo)V)']'''}  (18) 
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The  strengths  and  spaeing  of  the  pair  was  then  adjusted  sueh  that  a  stagnation 
streamline,  representative  of  the  model  volume,  was  formed.  The  stagnation 
streamline  is  loeated  at  the  x,  and  r  loeations  sueh  that  ¥=0. 


3.1.2  Wake  Simulation 

To  simulate  the  wake  of  the  propeller,  a  souree  or  sink  plaeed  at  the  propeller  loeation 
is  superimposed  with  uniform  flow  and  the  solid  body  souree  and  sink  pair.  The 
determination  of  whieh  singularity  sign  for  superimposition  with  uniform  flow  at  the 
propeller  location  depends  on  the  speed  of  the  wind  tunnel,  and  the  advance  ratio,  and 
thus  the  thrust  being  produced  by  the  propeller.  The  propeller  itself  produces  drag  for 
advance  ratios  above  a  critical  value  determined  by  the  pitch  of  the  propeller.  If  the 
propeller  is  unpowered,  the  advance  ratio  will  be  above  this  critical  value  and  thrust 
will  be  produced.  Therefore,  at  a  low  RPM  settings  and  low  wind  tunnel  speed,  the 
propeller  produces  drag,  rather  than  thrust,  and  the  wake  is  decelerated  when  the  net 
streamwise  force  on  the  propeller  and  nacelle  is  in  the  drag  direction.  In  this  case,  it 
is  appropriate  to  simulate  the  wake  as  a  source.  At  a  mid  power  setting  and  low  wind 
tunnel  speed,  the  wake  is  faster  than  the  surrounding  air.  A  sink  accurately  simulates 
the  wake.  At  a  high  power  setting,  the  wake  is  speeded  up  the  most  and  again  a  sink 
simulates  the  wake  accurately. 

By  the  momentum  theorem,  the  wake  will  be,  on  average,  faster  than  the  free  stream 
if  the  propeller  is  producing  a  net  thrust  and  the  wake  will,  on  average,  be  slower  than 
the  free  stream  if  the  propeller  is  producing  drag.  With  the  addition  of  a  source  or 
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sink  placed  at  the  propeller  loeation  for  wake  simulation  purposes,  a  singularity  of 
opposite  sign  and  equal  strength  must  be  placed  at  infinity  downstream. 

The  velocity  potential  for  a  point  souree  of  strength  Ows,  located  at  [x,y,z]  =  [-Xo,0,0] 
is 


<l)(x,y,z)  =  -  Ows/  {(47i)[(x+Xo)^+(y)^+(z)]'^^}  (19) 


The  veloeity  potential  for  a  point  sink  of  strength  Owk,  located  at  [x,y,z]  =  [-Xo,0,0]  is 


<l)(x,y,z)  =  Owk/  {(47i)[(x+Xo)^+(y)^+(z)^]'^^}  (20) 


The  corresponding  stream  function,  with  source  and  sink  locations  of  [x,r]  =  [-Xo,0] 
and  [x,r]  =  [xo,0],  respectively,  is 

¥(x,r)  =  -ows  (x+xo)  /  {2[(x-Xo)'  +  (21) 

'P(x,r)  =  Owk  (x-Xo)  /  {2[(x-Xo)^  +  (22) 


It  should  be  noted  that  the  stream  funetion  is  expressed  in  eylindrical  coordinates  due 
to  the  model  being  simulated  as  a  body  of  revolution  and  thus  being  axi-symmetrie. 
The  stream  function,  in  three-dimensions,  only  applies  to  axisymmetric  flow;  it  is  an 
approximate  description  of  axi-symmetric  body  in  a  rectangular  wind  tunnel. 
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3.2  Velocity  Increment  Calculation 


The  velocity  increment  induced  by  the  above  superimposition  of  uniform  flow,  solid 
body  source  and  sink  pair,  and  wake  source  and  sink  parr  was  calculated  next.  The 
velocity  induced  by  the  singularities  was  derived  from  the  velocity  potential  and 
stream  function  of  the  combination  of  the  singularities,  as  given  previously.  The 
induced  velocity  equations  due  to  the  singularities  follow. 


Induced  velocity  for  a  point  source  at  [x,y,z]  =  [xo,0,0] 

Au(x,y,z)  =  a(x-Xo)/  {(47i)[(x-Xo)^+(y)^+(z)^]^^^}  (23) 

Induced  velocity  for  a  point  sink  at  [x,y,z]  =  [xo,0,0] 

Au(x,y,z)  =  -  a(x-Xo)/  {(47i)[(x-Xo)^+(y)^+(z)^]^^^}  (24) 


Induced  velocity  for  uniform  flow 
Au(x,y,z)  =  Uoo 


(25) 


Since  the  model  was  mounted  in  the  center  of  the  tunnel,  all  singularities  lie  at  the 
center  of  the  tunnel.  Using  a  subscript  of  ss  to  indicate  a  solid  body  source,  a 
subscript  of  sk  to  indicate  a  solid  body  sink,  a  subscript  of  ws  to  indicate  a  wake 
blockage  source,  a  subscript  of  wk  to  indicate  a  wake  blockage  sink,  a  subscript  of  0 
to  indicate  the  coordinates  of  the  solid  body  singularities  ,  a  subscript  of  1  to  indicate 
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the  coordinates  of  the  wake  singularity  at  the  propeller,  and  a  subscript  of  2  to 
indicate  the  coordinates  of  the  wake  singularity  placed  far  downstream,  the  total 
induced  velocity  for  uniform  flow  superimposed  with  a  solid  body  source  and  sink 
pair,  a  point  source  located  at  the  at  the  propeller,  and  a  point  sink  far  downstream  is 


Au(x,y,z)  =  Uoo  +  ass(x+Xo)/  {(47i)[(x+Xo)^+(y)^+(z)^]^^^}  - 
ask(x-Xo)/  {(47i)[(x-Xo)^+(y)^+(z)^]^^^}  + 
aws(x-Xoi)/  {(47i)[(x-Xoi)^+(y)^+(zf]^^^}  - 

awk(x-Xo2)/  {(47i)[(x-Xo2)^+(y)^+(z)^]^^^}  (26) 


The  total  induced  velocity  for  uniform  flow  superimposed  with  a  solid  body  source 
and  sink  pair,  a  point  sink  located  at  the  at  the  propeller,  and  a  point  source  far 
downstream  is 


Au(x,y,z)  =  Uoo  +  ass(x+Xo)/  {(47i)[(x+Xo)^+(y)^+(z)^]^^^}  - 

ask(x-Xo)/  {(47i)[(x-Xo)^+(y)^+(z)^]^^^}  - 
awk(x-Xoi)/  {(47i)[(x-Xoi)^+(y)^+(z)^]^^^}  + 

aws(x-Xo2)  /  {(47i)[(x-Xo2)"+(y)'+(z)  ']'''}  (27) 
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In  this  test,  the  pressure  signature  was  taken  along  the  ceiling  of  the  wind  tunnel  at 
the  center  span.  The  coordinates  which  describe  the  pressure  taps  depend  both  on  the 
pressure  board’s  location  at  the  ceiling  as  well  as  the  location  of  the  pressure  taps 
along  the  wind  tunnel  test  section  length.  The  ceiling  location  and  this  “length” 
location  were  both  normalized  by  the  wind  tunnel  span.  This  results  in  the  first  tap  on 
the  pressure  board  being  located  at  [x/B,  y/B,  z/B]  =  [-0.6914,  0,  0.5].  The  last  tap  is 
located  at  [x/B,  y/B,  z/B]  =  [0.9336,  0,  0.5].  Pressure  taps  are  separated  by  x/B  = 
0.0125,  with  consistent  y/B  and  z/B  coordinates.  These  are  the  locations  at  which  the 
induced  velocity  must  be  calculated.  These  values  were  substituted  into  the  above 
appropriate  induced  velocity  equation  and  the  equation  was  further  divided  by  the 
given  uniform  flow  speed  and  square  of  the  tunnel  span  to  non-dimensionalize  the 
induced  velocity  equations.  The  velocity  increment  term,  Au  /  Uoo,  resulted. 

The  full  non-dimensional  equation  for  the  total  induced  velocity  for  uniform  flow 
superimposed  with  a  solid  body  source  and  sink  pair,  a  point  source  located  at  the  at 
the  propeller,  and  a  point  sink  far  downstream  is 


Au/  Uoo  =  1  +  ass(x/B+Xo/B)/  {(47i*  Uoo*B^)[(x/B+Xo/B)^  +  - 

ask(x/B-Xo/B)/  {(471*  Uoo*B^)[(x/B-Xo/B)^+  (z/B)^]^^^}  + 
aws(x/B-Xoi/B)/  {(471*  Uoo*B")[(x/B-Xoi/B)"+  (z/B)"]^^"}  - 
awk(x/B-Xo2/B) /  {(471*  Uoo*B")[(x/B-Xo2/B)"  +  (28) 
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The  full  non-dimensional  equation  for  the  total  indueed  veloeity  for  uniform  flow 
superimposed  with  a  solid  body  source  and  sink  pair,  a  point  sink  located  at  the  at  the 
propeller,  and  a  point  source  far  downstream  is 


Au/Uoo  =  1  +  ass(x/B+Xo/B)/  {(47i*  Uoo*B^)[(x/B+Xo/B)^+  (z/B)^]^^^}  - 

ask(x/B-Xo/B)/  {(471*  Uoo*b2)[(x/B-Xo/B)^+  (z/B)^]^^^}  - 
awk(x/B-Xoi/B)/  {(471*  Uoo*b2)[(x/B-Xoi/B)"  +  + 

aws(x/B-Xo2/B) /  {(471*  Uoo*b2)[(x/B-Xo2/B)'  +  (z/B) (29) 


To  completely  and  correctly  calculate  the  induced  velocity  due  to  the  wake  of  the 
propeller,  the  method  of  images  must  be  used  to  simulate  the  solid  boundaries  of  the 
wind  tunnel  walls.  This  is  presented  next. 


3,3  Method  of  Images 

The  method  of  images  is  a  technique  used  to  simulate  solid  boundaries,  such  as  wind 
tunnel  walls,  around  a  model  and  to  calculate  wall  interference  velocities  and 
blockage  effects.  In  this  case,  the  wake  of  the  propeller  is  simulated  with  either  a 
point  source  or  point  sink,  depending  on  the  working  state  of  the  propeller.  Using 
this  singularity  as  a  starting  point,  a  set  of  opposing  singularities  are  placed  at  an 
appropriate  distance  from  the  original  singularity  to  produce  a  zero  streamline  at  the 
location  of  the  solid  boundaries,  successfully  simulating  the  wind  tunnel  walls. 
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However,  these  additional  singularities  introduee  effeets  that  must  be  eaneelled  by 
plaeing  another  set  of  opposing  singularities  at  twice  the  distance  of  the  first.  These 
new  singularities  also  produce  additional  effects  and  another  set  of  opposing 
singularities  are  added  to  the  system.  This  obviously  becomes  an  infinite  process; 
point  sources  and  sinks  are  alternately  placed  at  appropriate  distances  Ifom  the 
original  singularity  to  simulate  the  solid  boundaries  and  while  producing  zero  net 
effect  due  to  these  singularities.  The  same  process  must  be  followed  to  simulate  wind 
tunnel  walls  at  the  singularity  of  opposite  strength  far  downstream.  In  theory,  sources 
and  sinks  should  be  placed  in  an  infinite  grid.  In  this  case,  22,500  singularities  were 
placed  around  the  original  singularity. 

Once  these  singularities  are  placed  and  the  image  system  is  established,  its  effect  on 
the  model  is  the  same  as  that  of  the  boundary  it  represents^.  Following  this 
establishment,  the  velocity  increment  produced  by  each  singularity,  including  those 
that  simulate  the  solid  body  and  wake,  is  calculated.  These  increments  are  summed, 
plotted  over  the  length  of  the  test  section,  and  compared  to  that  calculated  Ifom  the 
measured  pressure  signature.  These  increments  match  when  the  model  and  wake  has 
been  accurately  simulated  with  the  source  and  sink  singularities.  Since  the  model  is 
accurately  simulated  by  the  solid  body  source  and  sink  pair,  the  source  and  sink  pair 
that  simulate  the  wake  is  the  only  singularity  strength  that  is  adjusted  such  that  the 
velocity  increments  match. 
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3,4  Deduction  of  Singularity  Strength 


With  the  singularities  whieh  eorreetly  simulate  solid  bloekage,  wake  bloekage 
singularity  strength  and  sign  must  be  found  such  that  the  velocity  increment 
calculated  using  these  singularities  matches  the  velocity  increment  derived  from 
pressure  measurements  at  the  wall.  As  aforementioned,  the  sign  of  the  upstream 
singularity  is  easily  deduced  from  examining  the  velocity  increment  trend.  If  the 
increment  decreases  over  the  length  of  the  test  section,  the  propeller  is  simulated  by  a 
sink.  If  the  increment  increases  over  the  length  of  the  test  section,  the  propeller  is 
simulated  by  a  source.  As  expected,  the  magnitude  of  the  singularity  increases  with 
an  increase  in  magnitude  of  the  velocity  increment. 

For  each  data  run  that  was  selected  for  analysis,  singularity  strengths  were  found 
using  trial  and  error.  Initially,  the  mathematically  simulated  velocity  increments 
matched  well  in  shape  but  not  in  absolute  value.  As  is  described  in  Hackett,  et  al. 

On  the  majority  of  occasions,  the  test  section  will  be  too  short  in  relation  to 
the  model  and  its  wake  for  the  asymptote  to  be  well  established  at  both  ends. 
Designated  points  may  therefore  be  selected  at  the  upstream  and  downstream 
ends  of  the  signature.  On  the  first  attempt  to  match  the  curve,  these  points  are 
assumed  to  be  asymptotes.  Miss-distances  are  next  determined,  between  the 
designated  points  and  the  fitted  curve  and  the  curve  is  adjusted  so  that  the 
fitted  cure  intersects  the  designated  points.  This  correction  is  made  just  once 
and  is  not  iterated.^ 

Though  the  downstream  velocity  increments  appear  to  establish  asymptotic  behavior, 
the  upstream  asymptote  is  not  realized  and  determination  of  this  “miss  distance”  was 
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necessary.  Once  it  was  established  for  each  curve,  the  velocity  increment  calculated 
from  the  singularities  matched  the  experimentally  velocity  increments  very  well. 

For  purposes  of  demonstration,  data  will  be  presented  for  the  18”  x  8”  propeller 
experimentation.  All  other  data  may  be  found  in  Appendix  B.  Figures  2  through  9 
show  the  velocity  increment  calculated  using  the  singularities  to  simulate  solid  body 
blockage  and  wake  blockage  as  compared  to  experimental  data.  Each  figure  shows 
the  components  which  contribute  to  the  total  velocity  increment  at  the  ceiling  of  the 
tunnel,  however,  the  miss  distance  is  not  applied  to  the  components.  The  miss 
distance  is  applied  to  the  total  velocity  increment.  It  should  be  noted  that  the  velocity 
increment  caused  by  the  solid  body  blockage  was  virtually  zero  and  therefore  is  not 
shown.  During  the  1 8”  x  8”  propeller  testing,  there  was  no  high  speed,  high  RPM 
run  that  was  successfully  completed  at  NATF.  The  motor  overheated  on  each 
attempt;  there  are  eight  figures  below  showing  ah  other  speed  and  RPM  runs  selected 
for  analysis. 
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Velocity  Increment  vs  X/B  18x8  for  U=50ft/sec  and  Low  Power 
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Figure  2.  Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced 
velocity  increment  for  18”x8”  propeller  rotating  at  low  RPM  and  wind  tunnel  speed  of  50ft/s. 


Velocity  Increment  vs  X/B  18x8  for  U=50ft/sec  and  Medium  Power 
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Figure  3.  Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced 
velocity  increment  for  18”x8”  propeller  rotating  at  medium  RPM  and  wind  tunnel  speed  of  50ft/s. 


Velocity  Increment  vs  X/B  18x8  for  U=50fl/sec  and  High  Power 
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Figure  4.  Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced 
velocity  increment  for  18”x8”  propeller  rotating  at  high  RPM  and  wind  tunnel  speed  of  50ft/s. 


Velocity  Increment  vs  X/B  18x8  for  Ll=70ft/sec  and  Low  Power 
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Figure  5.  Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced 
velocity  increment  for  18”x8”  propeller  rotating  at  low  RPM  and  wind  tunnel  speed  of  70ft/s. 


Velocity  Increment  vs  X/B  18x8  for  U=70ft/sec  and  Medium  Power 
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Figure  6.  Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced 
velocity  increment  for  18”x8”  propeller  rotating  at  medium  RPM  and  wind  tunnel  speed  of  70ft/s. 


Velocity  Increment  vs  X/B  18x8  for  U=70fl/sec  and  High  Power 


juin/n  Biiep 


Figure  7.  Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced 
velocity  increment  for  18”x8”  propeller  rotating  at  high  RPM  and  wind  tunnel  speed  of  70fl/s. 


Velocity  Increment  vs  X/B  18x8  for  U=130ft/sec  and  Low  Power 
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Figure  8.  Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced 
velocity  increment  for  18”x8”  propeller  rotating  at  low  RPM  and  wind  tunnel  speed  ofl30ft/s. 


Velocity  Increment  vs  X/B  18x8  for  U=130fl/sec  and  Medum  Power 


Figure  9.  Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced 
velocity  increment  for  18”x8”  propeller  rotating  at  medium  RPM  and  wind  tunnel  speed  ofl30ft/s. 


Two  assessments  may  be  made  from  the  above  figures.  First,  the  methodology  used 
to  simulate  the  solid  body  and  the  wake  with  singularities  is  sound;  the  ealculated 
velocity  increment  matches  experimental  data  well  in  most  cases.  Some  cases  show 
data  scatter  which  cannot  be  identically  matched  by  the  mathematical  calculations. 
This  is  likely  a  result  of  data  scatter  inherent  in  the  data  collection  system.  Second, 
one  can  observe  that  as  the  power  setting  increases,  with  wind  tunnel  speed  held 
constant,  thrust  increases.  However,  as  wind  tunnel  speed  increases,  given  the  same 
power  setting,  thrust  produced  decreases.  This  is  an  expected  result. 

Figures  2  through  9  show  data  gathered  fromNATF.  Comparable  simulation  of  solid 
body  and  wake  blockage  were  not  performed  for  data  collected  from  GLMWT.  This 
is  because  there  was  no  blockage  effect  at  GLMWT.  This  can  be  seen  by  examining 
the  incremental  velocity  data  from  GLMWT  against  NATF.  Figures  10  through  12 
show  the  velocity  increments  from  GLMWT  as  compared  to  that  at  NATF.  The 
GLMWT  velocity  increments  show  that  there  is  no  appreciable  velocity  increment 
caused  by  the  propeller  thrust  at  any  wind  tunnel  velocity  or  RPM  setting.  There  is  a 
brief  drop  in  velocity  increment  at  x/B=0,  which  coincides  with  the  propeller  location. 
This  is  due  to  the  pressure  increase  through  the  propeller  disc.  Aside  from  this 
change,  the  velocity  increments  from  the  GLMWT  data  show  no  appreciable  slope. 

If  one  attempted  to  model  this  velocity  increment  with  singularities,  as  is  done  for  the 
NATF  velocity  increment,  one  would  discover  that  the  singularity  strength  would  be 
essentially  zero.  This  is  indicative  of  zero  blockage.  In  some  cases,  successful  runs 
at  NATF  were  not  able  to  be  successfully  completed  at  GLMWT  and  vice  versa. 
Corresponding  data  for  all  other  propellers  can  be  found  in  Appendix  B. 
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Additionally,  the  Glauert  correction  was  applied  to  the  GLMWT  test  data  to  prove 
that  no  blockage  corrections  were  necessary.  Figure  13-16  show  these  corrections. 

In  all  cases,  the  corrected  data  lies  directly  on  top  of  the  experimental  data,  or  very 
close  to  it.  In  one  case,  for  the  I8”x8”  propeller,  two  runs  were  completed.  These  are 
designated  by  the  diamonds  and  triangles.  The  repeatability  of  the  runs  was  not  as 
desired  and  caused  scatter.  This  is  shown  in  figure  14. 
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Velocity  Increment  vs  X/B  18x8  50ft/sec 
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Figure  10.  Velocity  increment  calculated  from  pressure  data  collected  at  NATF  and  GLMWT  for  the  18”x8”  pro 
tunnel  speed  of  50ft/s. 


Velocity  Increment  vs  X/B  18x8  70ft/s 
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Figure  1 1 .  Velocity  increment  calculated  from  pressure  data  collected  at  NATF  and  GLMWT  for  the  1 8”x8”  pro 
tunnel  speed  of  70ft/s. 


Velocity  Increment  vs  X/B  18x8  130ft/s 


Figure  12.  Velocity  increment  calculated  from  pressure  data  collected  at  NATF  and  GLMWT  for  the  18”x8”  pro 
tunnel  speed  of  70fl/s. 
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Ct  vs  J  for  18x8  GLMWT  &  NATF 
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Figure  14.  GLMWT  experimental  data  eompared  to  Glauert-eorrected  GLMWT  data. 


Ct  vs  J  for  22x8  GLMWT  &  NATF 
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3,5  Velocity  Increment  at  Wind  Tunnel  Centerline 


As  is  done  in  the  Haekett-Wilsden  method,  veloeity  inerements  at  the  tunnel 
eenterline  are  ealeulated.  This  is  done  simply  using  the  x/B  locations  of  the  pressure 
taps  while  setting  the  z/B  coordinate  at  which  the  velocity  increment  is  calculated  to 
the  centerline  location  of  z/B  =  0.  The  singularities  that  simulate  the  solid  blockage 
and  the  wake  blockage  are  removed  for  this  calculation,  and  only  the  increment 
caused  by  the  image  system  is  calculated.  The  increments  are  calculated  as  in 
equation  28,  for  each  singularity  in  the  infinity  of  images.  As  is  described  in  the 
method  of  images,  location  of  each  image  varies  in  y/B  and  z/B  coordinates. 
Additionally,  the  “miss  distance”  is  included  in  the  calculation  of  the  velocity 
increment  to  capture  the  aforementioned  effect  of  test  section  length. 

The  centerline  velocity  increments  are  shown  for  the  18”  x  8”  propeller  in  figures  17 
through  24.  Centerline  velocity  increments  for  all  other  propellers  can  be  found  in 
Appendix  B. 
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Figure  24.  Centerline  velocity  increment  for  18”  x  8”  propeller  at  13 


3,6  Calculation  of  Corrections 


Following  the  successful  modeling  of  the  model  and  propeller  wake,  the  classical 
Glauert  correction  was  calculated  as  described  in  chapter  1,  section  1.2.  As 
previously  stated,  this  calculation  provides  the  speed  at  which  the  experiment  should 
be  run  to  develop  the  forces  corresponding  to  those  that  would  be  obtained  in 
unrestricted  flow,  at  the  same  operating  conditions. 

It  should  be  noted  that  Glauert’s  correction  is  useful  both  for  positive  and  negative 
thrust  cases,  as  is  Sorensen  and  Mikkelsen’s.  However,  at  Cj  =  -1,  Glauert’s 
correction  becomes  singular.  Sorensen  and  Mikkelsen  state  that  the  correction 
formula  they  developed  allows  for  a  closed-form  solution  of  the  equations  of  axial 
momentum.  They  claim,  “With  the  new  approach,  there  is  no  longer  any  need  for 
introducing  approximate  correction  formulas  for  the  wind  speed  through  the  rotor 
disk.  Furthermore,  the  singularity  problem  is  solved  with  the  new  solution.^” 
Therefore,  corrections  were  calculated  only  both  positive  and  negative  thrust  runs  for 
comparison  between  methodologies.  The  Glauert  correction  calculated  from 
experimental  data  was  applied  to  negative  thrust  cases  and  plotted  against  the 
theoretical  corrections.  This  is  shown  below  in  figure  25.  Clearly,  the  Glauert 
corrections  obtained  with  experimental  data  correspond  with  theoretical  corrections  at 
identical  thrust  conditions. 
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Glauert  Corrections 
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Figure  25.  Glauert  corrections  plotted  against  experimental  data. 


To  complete  the  Haekett-Wilsden  eorreetion  method,  the  previously  ealeulated 
eenterline  veloeity  inerements  are  used.  The  value  of  the  centerline  veloeity 
inerement  is  added  to  the  datum  wind  tunnel  veloeity  to  determine  a  eorreeted 
veloeity,  as  is  done  in  both  the  Glauert  caleulation.  This  gives  the  effeetive  free 
stream  flow  needed  to  obtain  the  same  thrust  eoeffieients  in  free  air  as  those  obtained 
in  the  wind  tunnel. 

Lastly,  Sorensen  and  Mikkelsen  eorreetion  is  performed.  In  order  to  earry  out  this 
eorreetion,  the  singularities  and  their  strengths  whieh  sueeessfully  model  solid  and 
wake  bloekage  must  be  known.  This  is  aoeomplished  using  the  Haekett-Wilsden 
method.  Next,  to  ealeulate  the  intermediate  value,  u,  as  is  neeessary  in  solving  the 
first  equation  of  the  Sorensen  and  Mikkelsen  methodology,  the  eross  seetional  area  of 
the  slip  stream  far  downstream  must  be  ealeulated.  This  was  done  using  the 
eonservation  of  momentum.  Onee  the  eross  seetion  area  of  the  slip  stream  far 
downstream  is  known,  the  value  P  was  easily  found  and  used  in  equation  7,  along 
with  the  strength  of  the  propeller  singularity  to  find  u.  With  this  variable  ealeulated, 
the  five  momentum  equations  (10-14)  were  solved  and  all  other  variables  ealeulated. 
Finally,  the  Sorensen  and  Mikkelsen  eorreetion  formula,  equation  (15),  was 
performed.  Just  as  with  the  Glauert  and  Haekett-Wilsden  eorreetion,  this  formula 
provides  the  speed  at  which  the  experiment  should  be  run  to  develop  the  forees 
corresponding  to  those  that  would  be  obtained  in  unrestrieted  flow. 

Using  the  new  airspeeds  ealeulated  from  the  Glauert,  Haekett-Wilsden,  and  Sorensen 
and  Mikkelsen  eorreetion  formulas,eorreeted  advanee  ratios  were  ealeulated.  To 
eompare  the  correction  to  the  experimentally  obtained  data,  and  to  eompare  the  each 
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method  against  the  others,  experimental  thrust  eoeffieient  was  plotted  against  the 
experimental  advanee  ratio  and  corrected  points  were  super  imposed  on  this  plot. 
These  comparisons  are  shown  below  in  figures  26  through  29. 
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Ct  vs  J  for  14x8  GLMWT  &  NATF 
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Ct  vs  J  for  18x8  GLMWT  &  NATF 


Figure  27.  Data  points  corrected  for  wind  tunnel  blockage  plotted  against  experimental  data  for  the  18”x8”  propeller. 
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As  can  be  seen  in  the  24”  x  8”  propeller  data,  shown  in  figure  29,  there  appears  to  be 
two  distinet  thrust  coefficient  vs.  advanee  ratio  curves.  They  have  eaeh  been 
delineated  by  a  diamond  and  triangle  around  the  data  points.  The  eause  of  this 
“double  eurve”  has  been  explored  thoroughly.  Based  on  this  exploration,  the  most- 
likely  cause  of  the  double  curve  is  due  to  hysteresis  in  the  data  eollection  equipment 
and  scatter  in  the  data.  The  examination  of  possible  causes  will  be  addressed  later. 
Interestingly,  this  data  split  did  not  oceur  during  testing  of  any  other  propellers  or 
during  NATF  tests,  though  the  same  data  collection  technique  was  used. 


3.7  Pressure  Change  Predictions 

In  each  of  Glauert’s  and  Sorensen  and  Mikkelsen’s  derivation  of  their  respective 
correction  formula,  an  expression  is  derived  for  the  pressure  inerease  or  deerease  due 
to  the  total  foree  the  propeller  exerts  on  the  air  in  the  wind  tunnel  test  seetion.  In 
Glauert’s  derivation,  four  equations  and  four  unknowns  can  be  used  to  solve  for 
pressure  differential.  The  first  two  equations  are  equations  of  eontinuity: 


II 

(30) 

(C  -  S)  U2  =  CVo  -  Au 

(31) 

The  second  two  equations  are  axial  momentum  equations; 
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T  =  Sipui(ui  -  Vo)  +  (C  -  S)pU2(U2  -  Vo)  +  C(P2  -  Po) 


(32) 


Where  T  =  thrust 

Po  =  wind  tunnel  datum  pressure 
P2  =  pressure  far  downstream  of  propeller 
T  =  Sipu(ui  -  Vo)  +  (C  -  S)pu2(u2  -  Vo)  +  l/2Cp(Vo"  -  U2")  (33) 


Using  the  equations  of  eontinuity  to  eliminate  u  solve  for  ui  in  terms  of  U2  then 
substituting  into  equation  33,  the  pressure  differential  ean  be  subsequently  solved. 

Sorensen  and  Mtkkelsen  have  an  explieit  equation  for  predicting  the  pressure  increase 
or  decrease  caused  by  the  propeller.  It  is  obtained  by  applying  the  momentum 
theorem  on  the  wind  tunnel: 


T  +  (pi  -  Po)C  =  puiSi(ui  -  Vo)  -  pu2(C  -  Si)(Vo  -  U2) 


(34) 


All  quantities  are  known  and  the  pressure  differential  can  be  calculated  directly. 

These  pressure  differential  predictions  were  transformed  into  velocity  increment 
predictions  using  the  same  method  as  was  used  for  experimental  pressure 
measurement. 

Table  1  shows  side  by  side  comparison  of  the  calculated  pressure,  coefficient  of 
pressure,  and  velocity  increment  differentials  for  the  18”x8”  propeller.  Additionally, 
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figures  30  through  32  show  the  velocity  increment  differential  as  predicted  by 
Glauert’s  and  Sorensen  and  Mikkelsen’s  methodology. 
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Figure  30.  Predicted  incremental  velocity  differentials  plotted  against  velocity  increments  obtained  from  experimental  data  for  the 
18”x8”  propeller  and  wind  tunnel  velocity  of  50  fl/s. 
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Figure  31.  Predicted  incremental  velocity  differentials  plotted  against  velocity  increments  obtained  from  experimental  data  for  the 
18”x8”  propeller  and  wind  tunnel  velocity  of  70  ft/s. 
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Figure  32.  Predicted  incremental  velocity  differentials  plotted  against  velocity  increments  obtained  from  experimental  data  for  the 
18”x8”  propeller  and  wind  tunnel  velocity  of  130  ft/s. 


3,8  Discussion 


3,8,1  Corrected  NATF  Data 

As  stated,  the  purpose  of  this  thesis  was  to  collect  high  quality  powered  propeller  data 
to  which  select  correction  methodologies  may  be  applied.  Glauert  corrections, 
Hackett-Wilsden,  and  Sorensen  and  Mikkelsen  correction  methodologies  were 
explored  to  assess  the  validity  of  the  two  latter  against  the  widely  used  former 
correction  method  for  wind  tunnel  blockage. 

Due  to  blockage  effects  and  consistent  with  Glauert’s  correction,  higher  thrust  is 
expected  to  result  in  the  model  and  tunnel  configuration  with  higher  blockage. 
Glauert’s  correction  states  that  due  to  the  necessary  satisfaction  of  continuity,  the 
discontinuity  in  pressure  across  the  propeller  blade,  and  changes  in  air  velocity 
moving  outside  of  the  propeller  slipstream,  a  higher  thrust  is  developed  by  the 
propeller  in  the  wind  tunnel  than  otherwise  would  be  in  an  unrestricted  flow  of  the 
same  speed  with  the  same  propeller  rotation  rate  and  blade  angle.  This  is  true  when 
the  propeller  is  operating  such  that  positive  thrust  is  produced.  In  this  case,  the 
negligibly  small  blockage  ratio  at  GLMWT  results  in  essentially  unrestricted  flow 
around  the  propeller.  Conversely,  at  NATF,  the  blockage  is  high  and  more  thrust  is 
expected  to  be  developed  by  the  propeller  at  NATF  than  at  GLMWT,  at  the  same 
operating  conditions. 

Figures  26  through  29  provide  comparisons  for  the  correction  methodologies.  These 
figures  show  thrust  coefficient  vs.  advance  ratio  data  for  both  tunnels.  Superimposed 
on  this  data  is  the  corrected  NATF  data,  resulting  from  calculated  corrections  using 


77 


the  Glauert,  Hackett-Wilsden,  and  the  S0rensen  and  Mikkelsen  methodologies.  Three 
assessments  ean  be  made: 


1)  The  magnitude  of  the  corrections  increase  as  propeller  diameter  increases. 

2)  Agreement  between  Glauert,  Hackett-Wilsden,  and  Sorensen  and  Mikkelsen 
appear  to  be  good  at  low  positive  thrust  coefficients,  but  begin  to  diverge  at 
thrust  coefficients  of  about  1.2. 

3)  There  appears  to  be  better  agreement  between  Glauert  and  Hackett-Wilsden 
corrections  than  between  Glauert  and  Sorensen  and  Mikkelsen  or  Hackett- 
Wilsden  and  Sorensen  and  Mikkelsen. 

Data  from  each  propeller  will  be  addressed  separately.  Figure  26  shows  the 
coefficient  vs.  advance  ratio  data  for  the  14”x8”  propeller.  Clearly,  there  was  little 
blockage  effect  even  at  NATF,  as  the  area  ratio  was  only  6.7%.  The  corrected  data, 
indicated  by  the  black  outlined  points  on  Figure  26,  is  very  close  to  the  actual  data 
collected. 

The  data  collected  for  the  18”x8”  data,  shown  in  figure  27,  demonstrates  that  the 
thrust  measured  at  NATF  was  not  appreciably  greater  than  that  GLMWT,  as 
expected.  This  is  surprising  considering  the  blockage  for  the  18”x8”  propeller  is 
11%.  However,  this  does  not  imply  that  there  was  no  blockage  correction  needed. 
The  corrected  NATF  data,  shown  in  black  outlined  points,  does  correct  the  data. 
However,  it  brings  the  NATF  data  curve  below  that  of  the  GLMWT  data.  This  may 
lie  in  the  scatter,  but  it  does  show  the  beginning  of  an  unexpected  trend. 
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Even  though  the  bloekage  eorreetion  shows  unexpeeted  results,  among  the  three 
correetion  teehniques,  good  agreement  is  shown.  Again,  the  eorrespondence  of 
Glauert  and  Sorensen  and  Mikkelsen  eorreetions  diverge  at  a  thrust  eoeffieient  of 
about  1.2. 

It  should  be  noted  that  in  figure  27,  both  the  repeat  run  for  the  18”x8”  propeller  at 
GLMWT  is  shown  as  well  as  a  repeat  run  for  NATE.  The  NATE  repeat  run  is  shown 
by  pink  outlined  triangles.  The  repeatability  at  NATE  is  outstanding.  Additionally, 
the  Glauert-eorreeted  GEMWT  data  is  shown  for  fully  eorreeted  eomparison 
purposes. 

Eigure  28  shows  the  eoeffieient  of  thrust  vs.  advance  ratio  for  the  22”x8”  propeller. 
This  data  set  continues  to  show  unexpected  results.  Eirst,  the  thrust  coefficient  data 
IfomNATE  is  not  appreciably  greater  than  the  same  data  Ifom  GLMWT.  This  is 
unexpected,  as  blockage  effects  should  be  more  pronounced  with  the  large  blockage 
ratio  of  16.5%  that  the  larger  propeller  casues.  Additionally,  the  corrected  data 
greatly  reduces  the  thrust  coefficient  for  the  NATE  data.  This  is  a  continuation  of  the 
unexpected  trend  in  correction  results.  While  care  has  been  taken  in  examining  data 
reduction  and  calculation  techniques,  the  cause  of  the  unexpected  result  has  not  been 
found.  Possible  causes  will  be  addressed  in  a  following  section. 

Einally,  figure  29  shows  the  thrust  coefficient  vs.  advance  ratio  data  for  the  24”x8” 
propeller.  Again,  the  result  is  unexpected,  but  it  is  also  consistent  with  the  18”x8” 
and  22”x8”  data  correction  trend;  the  thrust  coefficient  data  Ifom  NATE  is  less  than 
that  at  GLMWT.  Additionally,  low  advance  ratio  data  was  impossible  to  collect  at 
NATE  as  the  motor  consistently  overheated  prior  to  a  successful  test  point  being 
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taken.  Unfortunately,  this  prohibits  eomparison  of  data  at  advance  ratios  less  than 
0.37.  The  NATF  data  appears  to  he  within  the  “lower  curve”  of  the  GLMWT  data  at 
advance  ratios  of  about  0.45  and  greater.  However,  the  trend  shows  that  it  begins  to 
diverge  below  this  advance  ratio.  Additionally,  as  in  the  18”x8”  and  22”x8”  propeller 
data,  the  corrected  NATF  data  further  reduces  the  thrust  coefficient  data  instead  of 
correcting  it  to  the  GLMWT  data. 


3.8.2  Pressure  Prediction  Data 

Table  1  shows  the  predicted  pressure  change  and  corresponding  coefficient  of 
pressure  and  velocity  increment  changes  for  the  18”x8”  propeller.  Additionally, 
figures  30  through  32  show  the  predicted  velocity  increment  change.  This  data  was 
calculated  using  intermediate  equations  in  the  Glauert  and  Sorensen  and  Mtkkelsen 
correction  methodologies. 

Each  of  these  figures  shows  a  poor  prediction  of  velocity  increment  change  using 
either  methodology.  There  appears  to  be  no  trend  in  terms  of  over  prediction  or  under 
prediction  of  the  velocity  increment  change.  Hacketf  s  suggestion  that  the  test  section 
length  may  cause  lack  of  asymptote  establishment  potentially  explains  one  possible 
contributor  to  these  incorrect  predictions.  However,  in  the  case  of  this  experiment,  it 
appears  that  the  downstream  asymptotic  value  is  certainly  established  and  the  test 
section  length  is  not  a  cause  of  the  inaccuracy  between  the  predicted  downstream 
pressure  values  and  the  actual  downstream  pressure  values.  This  is  an  area  for  further 
study. 
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3.8.3  Possible  Reasons  for  Unexpected  Data 


3.8.3. 1  “Double  Curve”  in  24”x8”  Propeller  Data 

The  double  curve  in  the  24”x8”  propeller  data  is  unexpected  and  is  not  duplicated  in 
any  other  test  run  at  either  tunnel.  Several  possibilities  for  this  double  curve  behavior 
have  been  explored; 

1)  Due  to  time  constraints,  an  RPM  counter  was  not  included  in  the 
wind  tunnel  model.  This  resulted  in  using  the  stroboscope  for 
measurement  of  RPM.  This  proved  to  be  very  time  consuming. 
Additionally,  it  was  possible  that  if  the  stroboscope  was  not  used 
properly,  the  wrong  RPM  could  be  measured.  Great  care  was  taken 
at  NATF  that  this  be  avoided.  However,  the  potential  for  this  to 
occur  was  not  realized  at  GLMWT  and  not  as  much  care  was  used 
in  measuring  RPM.  Therefore,  it  is  possible  that  the  RPM  value 
recorded  may  be  a  blade  passage  frequency  in  error.  However,  after 
all  data  reduction  and  plotting  was  accomplished,  a  study  was 
completed  to  assess  the  probability  of  RPM  error  being  the  source 
of  the  double  curve.  It  was  concluded  that  the  RPM  was  reliable 
and  correct  and  this  possibility  was  discarded. 

2)  Hysteresis  in  the  power  setting  mechanism  on  the  remote  control  is 
a  possible  cause  for  the  double  curve.  This  was  realized  as  a 
possibility  only  after  observing  the  unusual  behavior  and  there  is  no 
way  to  measure  the  probability  of  this  cause  as  the  remote  control 
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has  been  dismantled.  However,  it  very  well  may  be  the  eause  of  the 
error  and  is  assigned  medium  probability. 

3)  It  is  also  possible  that  this  double  curve  is  due  to  RPM  effects  of  the 
propeller.  This  was  not  explored  and  it  is  recommended  that  it  be 
an  area  of  further  study. 

4)  Tip  mach  number  effects  may  also  cause  this  double  curve. 
However,  this  has  been  assigned  low  probability  based  on  the  fact 
that  the  smaller  diameter  propellers,  operating  at  higher  RPMs  did 
not  display  this  behavior.  The  larger  diameter  propeller  would  have 
to  be  at  a  much  higher  RPM  to  cause  tip  mach  number  effects. 

5)  Due  to  the  results  shown  in  the  18”x8”  propeller  repeat  runs,  it  is 
possible  that  the  data  scatter  and/or  lack  of  high  quality  repeatability 
at  GLMWT  causes  data  anomalies.  A  quality  examination  of  the 
repeatability  was  not  performed  both  due  to  time  constraints  and 
due  to  the  lack  of  repeat  runs  to  provide  sufficient  data  for  this 
study. 


3, 8,3,2  Blockage  Corrections  causing  Unexpected  Results 

In  addition  to  the  double  curve  shown  in  the  24”x8”,  the  results  given  by  the  blockage 
corrections  at  NATF  are  unexpected.  The  result  desired  in  this  study  was  to 
successfully  correct  higher-thrust  NATF  data  to  GLMWT  data  using  the  three 
methodologies.  This  clearly  was  not  the  result.  It  is  difficult  to  understand  why  the 
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thrust  values  at  NATF  were  not  higher  than  those  at  GLMWT  as  is  expeeted.  Only  a 
few  possibilities  exist  for  this: 

1)  The  data  gathered  at  GLMWT  was  not  eorrect.  Based  on  present  and  past 
performance  at  GLMWT,  this  is  not  highly  probable.  GLMWT  is  a  well- 
known  and  well-respected  facility  that  continuously  provides  high  quality 
data  to  industry  and  research  groups.  The  balance  at  GLMWT  has  a  low 
uncertainty  of  only  ±0.015  lb,  in  axial  force,  and  continues  to  operate  well 
through  the  present  time. 

2)  The  data  gathered  at  NATF  was  not  correct.  This  has  a  bit  higher 
probability  than  at  GLMWT  based  on  the  higher  uncertainty  of  ±0.49  lb 
versus  only  ±0.015  lb  at  GLMWT.  This  approximate  half-pound,  if 
realized,  may  cause  enough  inaccuracy  in  the  data  to  cause  produce  the 
unexpected  results.  However,  as  at  GLMWT,  the  NATF  tunnel  is  widely- 
used  and  well-respected  for  high  quality  data.  It  is  unlikely  that  this  is  the 
cause  in  the  unexpected  data. 

Neither  of  the  above  cases  seem  likely  to  have  caused  the  unexpected  results.  It  is 
suggested  that  this  also  be  an  area  for  future  study. 
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CHAPTER  4  CONCLUSIONS 


The  method  of  mathematieally  simulating  the  solid  body  and  wake  bloekage 
suggested  by  Haekett-Wilsden  is  a  sound  method.  It  provides  for  further  calculation 
of  pressures,  using  other  methodologies,  and  velocity  increments  and  allows  one  to 
compare  the  expected  data  to  the  experimentally  obtained  data. 

For  the  14”x8”  propeller,  little  blockage  is  observed  in  the  NATF  tunnel.  However, 
the  blockage  corrections  increase  as  the  propeller  diameter  increases.  This  is 
expected.  However,  the  absolute  value  of  the  thrust  data  collected  at  NATF  is 
unexpectedly  low  for  all  propellers  from  the  18”x8”  propeller  and  larger 

The  Glauert  and  Sorensen  and  Mtkkelsen  correction  methodologies  appear  to  agree 
for  thrust  coefficients  up  to  about  1.2.  They  begin  to  diverge  above  this  value.  Since 
Sorensen  and  Mikkelsen  applied  their  correction  methodology  to  wind  turbines  and 
not  to  thrust-producing  propellers,  perhaps  there  is  an  unexpected  effect  when  the 
propeller  is  in  a  thrust-producing  state.  Because  the  Glauert  correction  has  been  used 
and  accepted  for  over  seventy- five  years,  this  author  believes  that  this  is  the 
correction  that  should  be  applied  over  the  Haekett-Wilsden  and  Sorensen  and 
Mikkelsen  correction  until  the  divergence  between  the  corrections  is  understood. 

The  Glauert  and  Haekett-Wilsden  corrections  appear  to  agree  up  to  a  thrust 
coefficient  of  approximately  2.  This  shows  better  agreement  than  between  Glauert 
and  Sorensen  and  Mikkelsen. 

It  is  interesting  to  note  that  although  the  corrections  for  the  18”x8”,  22”x8”,  and 
24”x8”  propellers  appear  to  be  incorrect,  all  three  methodologies  resulted  in 
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approximately  the  same  eorreetion.  This  eomparison  aeeomplishes  the  original  goal 
of  evaluation  of  two  new  eorreetion  methodologies  to  the  relatively  aeeepted  Glauert 
eorreetion. 


85 


CHAPTER  5  RECOMMENDATIONS  FOR  FUTURE  WORK 


It  is  recommended  that  in  future  propeller  data  collection  that  specific  attention  be 
paid  to  the  methodology  of  collection.  Specifically,  experimentation  can  be  executed 
to  examine  the  effects  of  varying  RPM  settings  with  the  wind  tunnel  speed 
consistently  set  at  a  low  or  high  value.  Avoidance  of  two  “curves”  in  the  data  is 
desired.  Exploration  of  possible  tip  mach  effects  and  RPM  effects  should  be 
conducted. 

Examination  into  the  cause  of  the  divergence  between  Glauert  and  Sorensen  and 
Mikkelsen  corrections  above  the  thrust  coefficient  value  of  1.2  should  also  be 
completed. 

Cause  of  the  unexpectedly  low  thrust  coefficient  values  for  the  18”x8”,  22”x8”,  and 
24”x8”  propeller  data  in  NATE  should  be  understood. 

Einally,  the  change  in  pressure  predictions  calculated  tfom  the  Glauert  derivation  and 
the  Sorensen  and  Mikkelsen  derivation  should  certainly  be  explored  further.  The  lack 
of  trends  and  accuracies  is  unexpected. 
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APPENDIX  A 


Scimitar  series  propellers  with  8”  geometric  pitch. 
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Propeller  balance  stand,  drill  guide  bushing,  and  drill  bit. 
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Dimensions  In  Inches  (millimeters) 


Rear  View 


Power  supply. 
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Remote  control  with  trim  pot  modification. 


Stroboscope  at  NATF. 
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Pressure  board  installed  at  NATF 


Sensing  equipment  for  pressure  tubes. 
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Motor  with  housing  and  ESC  (left)  and  as  installed  (right). 


Motor  and  housing  installed  on  support  structure  at  GLMWT. 
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Motor  and  housing  installed  on  support  structure  at  NATF. 
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APPENDIX  B 
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NATF  drag  tare,  motor  spinning  and  stationary. 


Velocity  Increment  vs  X/B  14x8  for  U=501t/sec  and  Low  Power 


odd  d  d  d 


CD 


jUin/n  eiiep 


Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  14”x8”  propeller  rotating  at  low  RPM  and  wind  tunnel  speed  of  50ft/s. 


Velocity  Increment  vs  X/B  14x8  for  U=50ft/sec  and  Medium  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  14”x8”  propeller  rotating  at  medium  RPM  and  wind  tunnel  speed  of  50ft/s. 


Velocity  Increment  vs  X/B  14x8  for  U=50ft/sec  and  High  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  14”x8”  propeller  rotating  at  high  RPM  and  wind  tunnel  speed  of  50ft/s. 


Velocity  Increment  vs  X/B  14x8  for  U=80ft/sec  and  Medium  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  14”x8”  propeller  rotating  at  medium  RPM  and  wind  tunnel  speed  of  80ft/s. 


Velocity  Increment  vs  X/B  14x8  for  U=80ft/sec  and  High  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  14”x8”  propeller  rotating  at  medium  RPM  and  wind  tunnel  speed  of  80ft/s. 


Velocity  Increment  vs  X/B  14x8  for  U=1101t/sec  and  Medium  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  14”x8”  propeller  rotating  at  medium  RPM  and  wind  tunnel  speed  of  1  lOft/s. 


Velocity  Increment  vs  X/B  14x8  for  U=11  Oft/sec  and  High  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  14”x8”  propeller  rotating  at  high  RPM  and  wind  tunnel  speed  of  1  lOft/s. 


Velocity  Increment  vs  X/B  22x8  for  U=50ft/sec  and  Low  Power 


fuin  /  n  eiiep 


Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  22”x8”  propeller  rotating  at  low  RPM  and  wind  tunnel  speed  of  50fl/s. 


Velocity  Increment  vs  X/B  22x8  for  U=50fl/sec  and  Medium  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  22”x8”  propeller  rotating  at  medium  RPM  and  wind  tunnel  speed  of  50fl/s. 


Velocity  Increment  vs  X/B  22x8  for  U=50ft/sec  and  High  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  22”x8”  propeller  rotating  at  high  RPM  and  wind  tunnel  speed  of  50ft/s. 


Velocity  Increment  vs  X/B  22x8  for  U=70ft/sec  and  Medium  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  22”x8”  propeller  rotating  at  medium  RPM  and  wind  tunnel  speed  of  70fl/s. 


Velocity  Increment  vs  X/B  22x8  for  U=70ft/sec  and  High  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  22”x8”  propeller  rotating  at  high  RPM  and  wind  tunnel  speed  of  70fl/s. 


Velocity  Increment  vs  X/B  22x8  for  U=90ft/sec  and  Medium  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  22”x8”  propeller  rotating  at  medium  RPM  and  wind  tunnel  speed  of  90fl/s. 


Velocity  Increment  vs  X/B  22x8  for  U=90ft/sec  and  High  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  22”x8”  propeller  rotating  at  high  RPM  and  wind  tunnel  speed  of  90ft/s. 


Velocity  Increment  vs  X/B  24x8  for  U=60ft/sec  and  Low  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  24”x8”  propeller  rotating  at  low  RPM  and  wind  tunnel  speed  of  60fl/s. 


Velocity  Increment  vs  X/B  24x8  for  U=60fl/sec  and  Medium  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  24”x8”  propeller  rotating  at  medium  RPM  and  wind  tunnel  speed  of  60fl/s. 


Velocity  Increment  vs  X/B  24x8  for  U=60ft/sec  and  High  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  24”x8”  propeller  rotating  at  high  RPM  and  wind  tunnel  speed  of  60ft/s. 


Velocity  Increment  vs  X/B  24x8  for  U=90ft/sec  and  Low  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  24”x8”  propeller  rotating  at  low  RPM  and  wind  tunnel  speed  of  90fl/s. 


Velocity  Increment  vs  X/B  24x8  for  L)=90ft/sec  and  Medium  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  24”x8”  propeller  rotating  at  medium  RPM  and  wind  tunnel  speed  of  90ft/s. 


Velocity  Increment  vs  X/B  24x8  for  U=120fl/sec  and  Low  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  24”x8”  propeller  rotating  at  low  RPM  and  wind  tunnel  speed  of  120ft/s. 


Velocity  Increment  vs  X/B  24x8  for  U=120ft/sec  and  Medium  Power 
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Velocity  increment  calculated  from  mathematical  simulation  model  and  propeller  plotted  against  experimentally  produced  velocity 
increment  for  24”x8”  propeller  rotating  at  medium  RPM  and  wind  tunnel  speed  of  120ft/s. 


Velocity  Increment  vs  X/B  14x8  50ft/s 
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Velocity  increment  calculated  from  pressure  data  collected  at  NATF  and  GLMWT  for  the  14”x8”  propeller  and  wind  tunnel  speed  of 
50ft/s. 


Velocity  Increment  vs  X/B  14x8  80ft/s 


n/nB}|9p 


Velocity  increment  calculated  from  pressure  data  collected  at  NATF  and  GLMWT  for  the  14”x8”  propeller  and  wind  tunnel  speed  of 
80ft/s. 


Velocity  Increment  vs  X/B  14x8  11  Oft/s 
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Velocity  increment  calculated  from  pressure  data  collected  at  NATF  and  GLMWT  for  the  14”x8”  propeller  and  wind  tunnel  speed  of 
11  Oft/s. 


Velocity  Increment  vs  X/B  22x8  50  ft/s 
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Velocity  increment  calculated  from  pressure  data  collected  at  NATF  and  GLMWT  for  the  22”x8”  propeller  and  wind  tunnel  speed  of 
50fl/s. 


Velocity  Increment  vs  X/B  22x8  70ft/s 
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Velocity  increment  calculated  from  pressure  data  collected  at  NATF  and  GLMWT  for  the  22”x8”  propeller  and  wind  tunnel  speed  of 
70ft/s. 


Velocity  Increment  vs  X/B  22x8  90ft/s 


Velocity  increment  calculated  from  pressure  data  collected  at  NATF  and  GLMWT  for  the  22”x8”  propeller  and  wind  tunnel  speed  of 
90ft/s. 


Velocity  Increment  vs  X/B  24x8  60ft/s 
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Velocity  increment  calculated  from  pressure  data  collected  at  NATF  and  GLMWT  for  the  24”x8”  propeller  and  wind  tunnel  speed  of 
60ft/s 


Velocity  Increment  vs  X/B  24x8  90ft/s 
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Velocity  increment  calculated  from  pressure  data  collected  at  NATF  and  GLMWT  for  the  24”x8”  propeller  and  wind  tunnel  speed  of 
90ft/s 


Glauert 

Deltau/U 

0.039374  I 

00 

00 

o 

o 

0.017195  1 

0.017957  1 

-0.0067  1 

-0.0304  1 

-0.0536 

Delta  Cp 

-0.080299 

00 

o 

o 

1 

LO 

00 

CD 

00 

o 

o 

1 

-0.036237 

00 

00 

CD 

CD 

CD 

O 

O 

00 

o 

o 

Delta  P 

-1.117251 

-0.595769 

-0.256562 

CD 

00 

CD 

00 

O 

CD 

1 

0.087988 

CD 

O 

CM 

CD 

O 

LO 

00 

o 

00 

o 

Sorensen 

Deltau/U 

-0.01207 

-0.00031 

-0.24999 

88699  0 

-0.00049 

-0.01371 

-0.09103 

Delta  Cp 

o 

o 

CN 

O 

O 

00 

CD 

O 

O 

O 

O 

0.43748 

-1.78851 

0.00098 

0.02723 

0.17378 

Delta  P  (psf) 

0.33393 

LO 

CD 

00 

O 

O 

O 

3.23609 

-5.10874 

0.00718 

0.07770 

0.49729 

O 

-0.49617 

-0.34072 

-0.27568 

-0.28819 

0.114477 

0.655731 

1 .905728 

Thrust 

(pounds) 

00 

00 

1 

-5.03 

-2.18 

-0.88 

0.9 

CM 

00 

00 

LO 

Density 

(psf) 

0.002287 

00 

00 

CM 

CM 

O 

O 

O 

0.002287 

0.002284 

00 

00 

CM 

CM 

o 

o 

o 

CD 

00 

CM 

CM 

O 

o 

o 

CM 

CD 

CM 

CM 

O 

o 

o 

Power  Setting 
(1=L,  2=M, 

3H) 

(M 

00 

CM 

oo 

CM 

00 

Velocity 

(ft/sec) 

110.3070 

109.8700 

(J> 

00 

CM 

CD 

00 

50.0125 

80.1783 

49.9617 

49.9712 

<D 

o 

c3 

<D 

zn 


c3 

<D 

c/5 

<D 

Q 

00 

4=) 

'T^ 

<D 

■4— » 

<!-) 

Oh 

C/) 

c3 

c/) 

<D 


ci 

<3^ 


<!-) 

O 

ci 


Q 

13 

> 

c3 


<D 

c/) 

zn 

<D 

<4-H  ;-4 

O  ^3^ 


r- 


Velocity  Increment  vs  X/B  14x8  50ft/s 
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incremental  velocity  differentials  plotted  against  veloeity  inerements  obtained  from  experimental  data 
and  wind  tunnel  velocity  of  50  ft/s. 


Velocity  Increment  vs  X/B  14x8  80ft/s 
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Predicted  incremental  velocity  differentials  plotted  against  velocity  increments  obtained  from  experimental  data  for  the  14”x8 
propeller  and  wind  tunnel  velocity  of  80  ft/s. 


Velocity  Increment  vs  X/B  14x8  llOft/s 
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Predicted  incremental  velocity  differentials  plotted  against  velocity  increments  obtained  from  experimental  data  for  the  14”x8 
propeller  and  wind  tunnel  velocity  of  1 10  ft/s. 
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Velocity  Increment  vs  X/B  22x8  50  ft/s 
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Predicted  incremental  velocity  differentials  plotted  against  velocity  increments  obtained  from  experimental  data  for  the  22”x8 
propeller  and  wind  tunnel  velocity  of  50  ft/s. 


Velocity  Increment  vs  X/B  22x8  70ft/s 
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Predicted  incremental  velocity  differentials  plotted  against  velocity  increments  obtained  from  experimental  data  for  the  22”x8 
propeller  and  wind  tunnel  velocity  of  70  ft/s. 


Velocity  Increment  vs  X/B  22x8  90ft/s 


Predicted  incremental  velocity  differentials  plotted  against  velocity  increments  obtained  from  experimental  data  for  the  22”x8 
propeller  and  wind  tunnel  velocity  of  90  ft/s. 
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Velocity  Increment  vs  X/B  24x8  60ft/s 
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Predicted  incremental  velocity  differentials  plotted  against  velocity  increments  obtained  from  experimental  data  for  the  24”x8 
propeller  and  wind  tunnel  velocity  of  60  ft/s. 


Velocity  Increment  vs  X/B  24x8  90ft/s 


Predicted  incremental  velocity  differentials  plotted  against  velocity  increments  obtained  from  experimental  data  for  the  24”x8 
propeller  and  wind  tunnel  velocity  of  90  ft/s. 


Velocity  Increment  vs  X/B  24x8  120ft/s 
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Predicted  incremental  velocity  differentials  plotted  against  velocity  increments  obtained  from  experimental  data  for  the  24”x8 
propeller  and  wind  tunnel  velocity  of  120  fl/s. 
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Intermediate  values  and  experimental  and  theoretical  Sorensen  and  Mikkelsen  corrections  for  the  14”x8”  propeller. 
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